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We squeeze unconditionally the collective spin of a dilute ensemble of laser-cooled 87 Rb atoms 
using their interaction with a driven optical resonator. The shape and size of the resulting spin uncer- 
tainty region are well described by a simple analytical model [M.H.S., I.D.L., V.V., arXiv:0911.3936 
through two orders of magnitude in the effective interaction strength, without free parameters. We 
deterministically generate states with up to 5.6(6) dB of metrologically relevant spin squeezing on 
the canonical 87 Rb hyperfine clock transition. 



Squeezed spin states [lrfl|j where a component of the 
total angular momentum of an ensemble of spins has 
less uncertainty 0, S| than is possible without quan- 
tum mechanical correlations attract interest for 
both fundamental and practical reasons. Fundamen- 
tally, they allow the study of many-body entanglement 
but retain a simple description in terms of a single col- 
lective angular- momentum variable 0, HI- Practically, 
they may be a means to overcome the projection noise 
limit on precision 0, 0, [HI, 14 1. Spin squeezing has 
been demonstrated using entanglement of ions via their 
shared motional modes [9|j , r epulsive interactions in a 
Bose-Einstein condensate [10|, or partial projection by 
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measurement 

In a companion paper [l5j we propose a cavity feed- 
back method for deterministic production of squeezed 
spin states using light-mediated interactions between dis- 
tant atoms in an optical resonator. This approach gener- 
ates spin dynamics similar to those of the one-axis twist- 
ing Hamiltonian H cx S z in Kitagawa and Ueda's original 
proposal [l| . Cavity squeezing scales to much higher par- 
ticle number than direct manipulation of ions i (but 
see Ref. [l6[ for a potentially scalable approach) and 
employs dilute ensembles rather than dense condensates 
of interac ting atoms (lOj . Unlike measurement- based 
squeezing [ 1 lL Hij]. it unconditionally produces a known 
squeezed state independent of detector performance. 

Here we implement cavity squeezing for the canonical 
\F = l,mp = 0) o \F = 2, nip = 0) hyperfine clock 
transition in 87 Rb atoms, achieving a 5.6(6) dB improve- 
ment in signal-to-noise ratio 0, [J. To our knowledge, 
this is the largest such improvement to date. More- 
over, the shape and orientation of the uncertainty re- 
gions we observe agree with a straightforward analytical 
model [15| . without free parameters, over two orders of 
magnitude in effective interaction strength. 

Our scheme, similar in spirit to the proposal of 
Ref. [TtJ , relies on the repeated interaction of the atomic 
ensemble with light circulating in an optical resonator, as 
illustrated in Fig.[TJ We label the two relevant eigenstates 
(clock states) of each one of No atoms as the spin-up and 
spin-down states of a spin- 1/2 Sj, and define a total spin 
S = 2~2i s i- Its z component corresponds to the popu- 
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FIG. 1: Cavity squeezing [H^ : (a) The atoms are trapped in 
a standing-wave dipole trap inside an optical resonator, (b) 
The probe laser is detuned from cavity resonance by half a 
linewidth, so that atom-induced shifts of the cavity frequency 
change the transmitted power, (c) The cavity is tuned halfway 
between the optical transition frequencies for the two clock 
states, (d) The ^-dependent light shift shears the circular 
uncertainty region of the initial coherent spin state (red cir- 
cle) into an ellipse (dotted). Photon shot noise causes phase 
broadening that increases the ellipse area (solid). The illus- 
tration is for a modest shearing Q = 3 (see text). 



lation difference between clock states and its azimuthal 
angle corresponds to their relative phase. For a given 
total spin magnitude S = \S\ < Sq = Nq/2 and a given 
permutation symmetry of the ensemble, the set of possi- 
ble collective states forms a Bloch sphere. 

The coupling of the atoms to the resonator manifests 
itself both as a differential light shift of the clock states 
which causes the Sj to precess about the z axis, and as a 
modified index of refraction which shifts the cavity reso- 
nance frequency. If a resonator mode is tuned halfway be- 
tween the optical transition frequencies for the two clock 
states [Fig. [He)], the atomic index of refraction produces 
opposite frequency shifts of the mode for atoms in each of 
the states, yielding a net shift Aw r /« = 4>iS z /2 <C 1 pro- 
portional to the population difference 2S Z . Here n is the 
linewidth of the resonator and <pi is the spin precession 
angle per photon transmitted through the resonator. The 
resonator is driven by a probe laser with fixed incident 
power at a detuning k/2 so that this mode frequency 
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shift changes the average number of photons transmit- 
ted by Ap = po4>iS z from its value po in the absence of 
atoms. As the intracavity power is ^-dependent, so is 
the light shift, which produces a precession of each spin 
through an angle <j>(S z ) = QS z /So. The state of each 
atom now depends, through S z , on that of all other atoms 
in the ensemble. The shearing strength Q = S$po4>\ is a 
dimensionless measure of the light-mediated interaction 
strength. In particular, a coherent spin state prepared 
on the equator of the Bloch sphere (an uncorelated state 
with (S) = Sx and AS 2 = AS 2 = 5/2) has its circular 
uncertainty region sheared into an ellipse with a short- 
ened minor axis [Fig.QJd)] fl5l |. 

Two fundamental decoherence mechanisms counteract 
the unitary evolution which squeezes the spin uncer- 
tainty. The first is photon shot noise: the intracav- 
ity light field, driven by a coherent input and decaying 
via the cavity mirrors, is not in a photon number state 
and produces an uncertain light shift. This uncertainty 
leads to irreversible phase broadening A(j> 2 = pQ(j>\/2 = 
Q/(2Sq). The squeezed variance, which would lessen as 
Q~ 2 if the dynamics preserved the area of the uncertainty 
region, therefore only decreases as Q~ x flit . 

The second decoherence process is photon scattering 
into free space. Scattered photons that reveal the state 
of individual atoms spoil the ensemble's coherence, while 
Raman scattering, which changes the atoms' internal 
state at random, increases the spin variance (l8| . In our 
system, Rayleigh scattering occurs at the same rate for 
the two clock states, does not reveal the atomic state, 
and so does not harm the coherence [l9j. At most 2.3% 
of the atoms undergo Raman scattering for o ur p aram- 
eters, causing added noise and decoherence [18( much 
smaller than those from technical sources. 

Finally, the coherent shearing action ceases to reduce 
the minimum spin variance once the uncertainty region 
becomes elongated enough that the curved geometry of 
the Bloch sphere becomes important [l|. Such curvature 
effects are visible in our data for large values of the shear- 
ing strength Q. 

For the experimental demonstration, up to A^ot = 
5 x 10 4 atoms of 87 Rb (with excited-state decay rate 
r = 2ir x 6.065 MHz) are confined in a standing-wave 
optical dipole trap inside a Fabry-Perot resonator of 
linewidth k = 2tt x 1.01(3) MHz. Details of the ap- 
paratus are given in Ref. [11|. The atoms are coupled 
to the resonator with a position-dependent dimension- 
less cooperativity r/(r) = 4o(r) 2 / kT, where 2g(r) is the 
vacuum Rabi frequency [20j. We define an effective co- 
operativity r\ and atom number Nq for an equivalent 
uniformly-coupled system so that the spin variance of 
a coherent spin state, measured via the resonator fre- 
quency shift, is AS 2 = S /2 = N /4. The effective 
quantities must satisfy NqT] = N to t(r](r)) e , where () e de- 
notes an average over the ensemble, in order to repro- 
duce the observed average frequency shift. Reproducing 
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FIG. 2: Normalized variance a 2 as a function of rotation an- 
gle a about the mean spin direction for states prepared with 
shearing Q = (black circles), Q = 1.2 (red squares), Q — 7.7 
(green triangles) and Q — 30.7 (blue diamonds). The curves 
are cosine fits. Statistical error bars are comparable to the 
symbol size. The shapes of the corresponding uncertainty re- 
gions are illustrated below the plot. Inset: Observed variance 
AS 2 of the initial state as a function of So- The line is the 
projection noise limit as determined from cavity parameters. 



the projection-noise-induced variance of the cavity shift 
requires NqT] 2 = N to t(r](r) 2 } e . These two constraints 
impose the definitions rj = (rj(r) 2 ) e / '(rj(r)} e = 0.139(5) 
and TVo = Ntot(v}e/v- The effective total spin has 
S z = (N2—Ni)/2, where Ni^ are the analogously-defined 
effective populations of the clock states. 

We prepare an initial coherent spin state by optically 
pumping the atoms into \F = l,mp = 0} ((S) = —Sqz) 
and applying a microwave tt/2 pulse to yield {S} = SqX. 
The squeezing is performed by two pulses of 780 nm light 
detuned 3.18(1) GHz to the blue of the \5 2 S 1/2 ,F = 2) f> 
|5 2 P 3 / 2 , F' = 3} transition and 500 kHz to the blue of the 
cavity resonance. This yields a single-photon phase shift 
<t>i = {2/2,)rjT/8 = 171(6) prad, where 5 = 2ir x 3.29 GHz 
is the effective detuning from the 5 2 P 3/ / 2 manifold with 
oscillator strength 2/3. Between the two optical pulses, 
each of which lasts 50 ps ^> and contains up to ~ 10 5 
photons, is a composite (SCROFULOUS [2l[) microwave 
7r pulse, forming a spin echo sequence. The spin echo can- 
cels the spatially inhomogeneous phase shift caused by 
the po photons transmitted on average through the res- 
onator but preserves the shearing effect. We measure S z 
using this same sequence but with stronger optical pulses, 
with 10 6 photons transmitted on average. The transmit- 
ted fraction of these pulses, measured on an avalanche 
photodiode, reveals the cavity resonance frequency shift 
and hence S z [llT |. 

We observe the shearing by rotating the state through 
an angle a with a microwave pulse about the axis of its 
mean spin vector and recording the variance AS" 2 of a 
subsequent measurement of S z \ a over a series of 100 iden- 
tical preparations. The measured variance is normalized 
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FIG. 3: o"m ax (top panel) is the normalized maximum variance, 
°min (middle panel) the normalized minimum variance, and Qo 
(bottom panel) the rotation angle for minimum variance for 
each measured ellipse as a function of shearing strength Q. 
Statistical error bars are given for o^im an d are smaller than 
the symbols for and ag. The blue dashed curves are 

theoretical predictions for an ideal system, while the black 
lines are predictions including separately measured technical 
imperfections. Shearing was varied by adjusting the photon 
number po « 2200 x Q. 



to projection noise cr 2 (a) = 2AS^/So where So is de- 
termined from our atom number calibration, based on 
a first-principles calculation using accurately measured 
cavity parameters 

[PH . The inset to Fig. [2] shows the ob- 
served A S 2 of the initial spin state as a function of So, 
as well as the calculated projection noise limit. 

Typical data of AS^ are displayed in Fig. [5J As the 
state is rotated the variance dips below projection noise 
as the minor axis of the ellipse is aligned with z and 
then increases beyond it as the major axis in turn rotates 
towards z. The variation of AS^ with angle is sinusoidal 
with a period 7T, as it must be for any distribution of 
S y -S z fluctuations. We record such curves over a range 
of photon numbers, corresponding to increasing shearing 
strength Q, keeping the effective atom number constant 
at 2So 3.2 x 10 4 . We compare the fitted phase and 
minimum and maximum variance of each sinusoid to the 
predictions of our model , briefly described below. 

Neglecting scattering, the initial S z distribution with 
AS Z = Sq/2 is unaffected by the shearing, while the S y 
variance is modified 1 1 511 
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and an S y -S z correlation is introduced 



W = (S y S z + S Z S V ) = tiQSc 



(1) 



(2) 



The expressions given here are approximations to those of 
the companion theory paper |15j . with which they agree 
to within 0.01% for our parameters. They include ad- 
ditional correction factors £, 7, and S/Sq to account for 
technical imperfections. £ = (d£/dcjp)«;/(2£) is the log- 
arithmic derivative of the Lorentzian resonator transmis- 
sion £ with respect to fractional probe detuning lu p /k. 
Ideally, uj p = k/2 and £ = 1 exactly. As we do not main- 
tain precisely this detuning, 0.97 < £ < 1 for the data 
presented here. The variance of the intracavity probe 
power in the absence of atoms, expressed as a multiple of 
photon shot noise, is 7 = l+2Ap 2 /p = l+po/(8x 10 4 ) = 
1 + Q/37, where Apj is the additional variance in the 
transmitted photon number caused by independently de- 
termined fractional light noise. Finally, dephasing re- 
duces the effective radius of the Bloch sphere from So to 
S, measured from the envelope amplitude of a Rabi nu- 
tation curve. For modest shearing Q < 20 we maintain 
a Bloch sphere radius S > O.8OS0, but when Q reaches 
200 the radius is reduced to S = 0.48S . 

Equations [TJ and [2] together yield a prediction for the 
observed normalized variance: 
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A cos (2a - 2a )] + 



(3) 



where V± = AS 2 ±AS 2 , A 



V 2 + W 2 , and the ro- 
tation angle which minimizes the variance is tan(2ao) = 
W/V-. The additional variance of our imperfect read- 
out <j 2 q — 0.13 is determined by comparing successive 
measurements of the same state as in Ref. [ll| . 

Figure [3] shows the predictions of maximum and mini- 
mum variance dmax = (V+ ± A) /So + o 2 Q and orientation 

min 

ao as black lines, together with the data points extracted 
from the cosine fits. It is remarkable that a simple an- 
alytical model, without free parameters, provides good 
predictions of the shape, area and orientation of the un- 
certainty region for values of the shearing strength Q that 
span a factor of 200. Note the increase of the minimum 
variance cr min for large shearing as the curvature of the 
Bloch sphere becomes important. For comparison, Fig. [3] 
also includes the model predictions without technical cor- 
rections (S = So, 7 = £ = 1, <J 2 = 0) as blue dashed 
curves. The maximal noise cr 2 ^^ and the ellipse angle a > 
dominated by the shearing-induced broadening along y, 
are insensitive to technical effects, but the minimum vari- 
ance cr mm is strongly affected by readout noise a 2 Q . This 
noise is due to a combination of finite quantum efficiency 
and avalanche noise of the photodetector together with 
Raman scattering which limits the number of photons 
in the readout measurement. It could be suppressed by 
using a different photodetector to remove the avalanche 
noise and by performing the readout near-detuned to a 
cycling transition to suppress Raman scattering [18|. 

To determine whether the reduced spin noise cr 2 ain al- 
lows a gain in spectroscopic precision, we must establish 
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FIG. 4: Metrologically-relevant squeezing £ (red solid circles) 
and signal contrast C (green diamonds, inset) as a function 
of shearing strength Q. 



a signal-to-noise ratio by comparing it to the mean spin 
signal \(S)\ @|. \(S)\ is reduced below Sq by dephas- 
ing, which shortens the Bloch vector, and by shot-to- 
shot phase fluctuations which reduce the average projec- 
tion of S along its mean direction. We specify the sig- 
nal strength by a contrast C = \(S)\/Sq, measured from 
the mean amplitude of Rabi oscillations with the sheared 
state as input, and plotted in the inset to Fig. 01 The sig- 
nal contrast in the absence of squeezing light C ln = 0.80 
is limited by dephasing from the lock light used to stabi- 
lize the cavity length 

Figure [4] shows the metrological squeezing parame- 
ter C = 2|<S)| in /(KS)| 2 /AS2) = a^C-JC 2 , which 
compares the squared signal-to-noise ratio for an ideal 
projection-noise-limited measurement using the initial 
spin signal | ((S) |i n = C ln So to that for the minimum ob- 
served variance and corresponding signal \ (S)\. Q < 1 in- 
dicates an improvement in signal-to-noise ratio unattain- 
able without entanglement [2|. For So = 1-6 x 10 4 and 
Po = 4.1 x 10 4 (Q = 19), we reduce the spin noise by 
a factor c^in = —6.7(6) dB. At this photon number 
our contrast is C = 0.78(2), so that we demonstrate a 
£ = 5.6(6) dB improvement in potential measurement 
precision over that of the initial uncorrelated state. 

By subtracting our independently measured readout 
noise ofo fr° m a mim we infer that states prepared by a 



shearing Q = 19 have an intrinsic spin variance that is a 
full 10(1) dB below the projection noise limit. However, 
it is the observed and not the intrinsic variance that de- 
termines the precision of a spectroscopic measurement, 
and we use the former in calculating the spin noise re- 
duction and the squeezing £. 

In conclusion, we have demonstrated a method for de- 
terministically generating squeezed states using switch- 
able light-mediated interactions in a dilute ensemble of 
otherwise non-interacting atoms. Our model predicts the 
size and shape of the uncertainty region when technical 
effects are included. We hope to observe states with sub- 



stantially lower spin noise by improving our readout, and 
to demonstrate their use in an atomic clock. 

This work was supported in part by the NSF, DARPA, 
and the NSF Center for Ultracold Atoms. M. H. S. ac- 
knowledges support from the Hertz Foundation and NSF. 
I. D. L. acknowledges support from NSERC. 



[1 
[2 

[3 

[4 

[5 

[6 

[7 

[8 

[9 

[10 
[11 

[12 

[13 
[14 

[15 

[16 
[17 

[18 
[19 



[20 
[21 



M. Kitagawa and M. Ueda, Phys. Rev. A 47, 5138 (1993). 

D. J. Wineland, J. J. Bollinger, W. M. Itano, F. L. Moore, 

and D. J. Heinzen, Phys. Rev. A 46, R6797 (1992). 

D. J. Wineland, J. J. Bollinger, W. M. Itano, and D. J. 

Heinzen, Phys. Rev. A 50, R67 (1994). 

A. S. S0rensen and K. M0lmer, Phys. Rev. Lett. 86, 4431 

(2001). 

A. S0rensen, L.-M. Duan, J. I. Cirac, and P. Zoller, Na- 
ture 409, 63 (2001). 

A. Kuzmich, N. P. Bigelow, and L. Mandel, Europhys. 
Lett. 42, 481 (1998). 

A. Kuzmich, L. Mandel, and N. P. Bigelow, Phys. Rev. 
Lett. 85, 1594 (2000). 

T. Takano, M. Fuyama, R. Namiki, and Y. Takahashi, 
Phys. Rev. Lett. 102, 033601 (2009). 
V. Meyer, M. A. Rowe, D. Kielpinski, C. A. Sackett, 
W. M. Itano, C. Monroe, and D. J. Wineland, Phys. 
Rev. Lett. 86, 5870 (2001). 

J. Esteve, C. Gross, A. Weller, S. Giovanazzi, and M. K. 

Oberthaler, Nature 455, 1216 (2008). 

M. H. Schleier-Smith, I. D. Leroux, and V. Vuletic, 

Reduced- quantum-uncertainty states for an atomic clock 

(2008), arXiv:0810.2582, to be published in Phys. Rev. 

Lett.. 

J. Appel, P. J. Windpassinger, D. Oblak, U. B. Hoff, 
N. Kjasrgaard, and E. S. Polzik, Proceedings of the Na- 
tional Academy of Sciences 106, 10960 (2009). 
A. Andre, A. S. S0rensen, and M. D. Lukin, Phys. Rev. 
Lett. 92, 230801 (2004). 

G. Santarelli, P. Laurent, P. Lemonde, A. Clairon, A. G. 
Mann, S. Chang, A. N. Luiten, and C. Salomon, Phys. 
Rev. Lett. 82, 4619 (1999). 

M. H. Schleier-Smith, I. D. Leroux, and V. Vuletic, 
Squeezing the collective spin of a dilute atomic ensemble 
by cavity feedback (2009), arXiv:0911.3936, to be pub- 
lished in Phys. Rev. A. 

D. Leibfried, M. D. Barrett, T. Schaetz, J. Britton, 
J. Chiaverini, W. M. Itano, J. D. Jost, C. Langer, and 
D. J. Wineland, Science 304, 1476 (2004). 
M. Takeuchi, S. Ichihara, T. Takano, M. Kumakura, 
T. Yabuzaki, and Y. Takahashi, Phys. Rev. Lett. 94, 
023003 (2005). 

M. Saffman, D. Oblak, J. Appel, and E. S. Polzik, Phys. 
Rev. A 79, 023831 (2009). 

R. Ozeri, C. Langer, J. D. Jost, B. DeMarco, A. Ben- 
Kish, B. R. Blakestad, J. Britton, J. Chiaverini, W. M. 
Itano, D. B. Hume, et al., Phys. Rev. Lett. 95, 030403 
(2005). 

H. J. Kimble, Physica Scripta T76, 127 (1998). 

H. K. Cummins, G. Llewellyn, and J. A. Jones, Phys. 
Rev. A 67, 042308 (2003). 



